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ABSTRACT 

Context. The first detection of a gamma ray line with an energy of about 500 keV from the center of our Galaxy dates back to the early 
seventies. Thanks to the astrophysical application of high spectral resolution detectors, it was soon clear that this radiation was due 
to the 511 keV photons generated by electron-positron annihilation. Even though the physical process are known, the astrophysical 
origin of this radiation is still a mystery. 

Aims. The spectrometer SPI aboard the INTEGRAL gamma-ray satellite has been used to produce the first all-sky map in light of the 
511 keV annihilation, but no direct evidence of any 511 keV galactic compact objects has been found. Owing to its moderate angular 
resolution, these SPI data are still compatible with a distribution of point sources clustered in the bulge of our Galaxy. Thanks to the 
fine angular resolution and the large field of view the IBIS imager on the INTEGRAL satellite gives us the unique opportunity to 
search for a possible 511 keV line from point sources associated to known objects, such as X-ray binaries, or supernovae, or even new 
ones. 

Methods. We present the first deep IBIS 511 keV all sky map, obtained by applying standard analysis to about 5 years of data. Possible 
511 keV signals are also searched over hour-day-month timescales. The IBIS sensitivity at 51 1 keV depends on the detector quantum 
efficiency at this energy and on the background. Both these quantities were estimated in this work. 

Results. We find no evidence of Galactic 511 keV point sources. With an exposure of 10 Ms, in the center of the Galaxy, we estimate 
a 1.6 X \Q^* phcnT- ^ ' flux 2 sigma upper limit. Similar limit is given in a wide area in the Galactic center region with similar 
exposures. The IBIS 511 keV flux upper limits for microquasars and supernova remnants detected in the hard X domain {E > 20 keV) 
are also reported. 

Conclusions. Our results are consistent with a diffuse e*e^ annihilation scenario. If positrons are generated in compact objects, we 
expect that a significant fraction of them propagate in the interstellar medium before they are annihilated away from their birth places. 

Key words. INTEGRAL - data analysis - e-l-e- annihilation 



1. Introduction 

Positrons can be observed through the detection of gamma rays 
produced by the electron-positron annihilation. This process can 
proceed in two ways, by direct annihilation or by the formation 
of positronium, i.e. bound states consisting of an electron and 
a positron. Twenty-five percent of the positronium atoms have 
antiparallel spins (para-positronium, p-Ps) and 75 % have paral- 
lel spins (oitho-positronium, o-Ps). It is known that p-Ps has a 
lifetime of 1.25 x 10"'" s and decays into two gamma rays each 
with 511 keV energy, while o-Ps has a lifetime of 1.5 x 10"^ s 
and three gamma rays are emitted, the maximum energy being 
511 keV. By evaluating the ratio between the line and the contin- 
uum component in the spectrum, one can estimate the fraction of 
positrons that decay by formation of positronium. 

The discovery of positrons in our Galaxy dates back to the 
seventies, when in a balloon experiment 'Johnson et al. 



1972) 1 



detected gamma ray line emission at the energy of 476 + 26 keV 
from the Galactic center. However, due to the low energy resolu- 
tion of the detectors, the physical origin of this emission was un- 
clear. A few years later, thanks to the advent of high-resolution 
spectrometers, this radiation was identified as a 511 keV nar- 
row line (FWHM < 3 keV) produced by electron-positron an- 
nihilation ( Leventhal et al.|[l978 i. Observations of the galactic 
annihilation radiation continued progressively until today with 



the new observations made with the INTErnational Gamma- 
Ray Astrophysics Laboratory (INTEGRAL). The INTEGRAL 
satellite, launched from Baikonur in Kazakhstan on October 17, 
2002, is an ESA mission in cooperation with Russia and the 
United States. 

The 511 keV flux (~ 10^^ phcm^^ s^^) measured by the 
spectrometer SPI aboard the INTEGRAL satellite shows that 
~ lO'*^ positrons per second annihilate in the bulge of our 
Galaxy ( Weidenspointner et al.|200^b] l. The origin of this large 
number of positrons is unclear. Indeed, although the physical 
processes responsible for producing positrons are quite well un- 
derstood, it is hard to find a picture that explains both this high 
rate of annihilation and the positron distribution in the Galaxy. 
The important observational constraints obtained so far can be 
divided into three main categories: (1) the spectral 511 keV line 
properties that give information on the environment where the 
positrons annihilate; (2) the spatial distribution of the emission 
(see e.g. [Weidenspointner et al. 2008a), which must be corre- 
lated with the possible positron sources candidates; and (3) the 
continuum galactic soft y-ray emission, which constrains the 
positron injection energy. Concerning points (1) and (2), the 
SPI data show that the 511 keV line width is on the order of 
keV, and the emission, contrary to what happens at other wave- 
lengths, is mainly coming from the Galactic bulge, in a region 
spatially extended on a scale of about 8 degrees. The measured 
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oitho-positronium continuum flux yields a fraction of positron- 
ium of (96.7 ± 2.2)% (Jean et al. 2006 1. Assuming a Galactic 
center distance of 8.5 kpc, the bulge size corresponds to approx- 
imately 1.2 kpc. Concerning point (3), the lack of detection of a 
y-ray continuum by COMPTEL/EGRET due to in-flight annihi- 
lation puts a severe constraint on the positrons injection energy. 
A calculation of the y-ray spectrum produced by the positron 
annihilation while they are still relativistic has shown that the 
emerging photon flux would exceed the COMPTEL/EGRET 
flux limit unless the positrons are injected into the interstellar 
medium with an energy below ~ 3 MeV (Beacom & Yiiksel 
[2006). An additional observational constraint is the flux upper 
limit on possible point sources. Possibly due to positron an- 
nihilation, transient broad-line emission has been reported in 
a few cases. High-energy emissions characterized by a broad 



line i n the 300-600 keV range from IE 1740.7-2942 (|Sunyaev 
et al . 1991; Bouc het et al.||1991 ) and Nova Muscae ( |Gilfanov 



et al.^199 1) have been detected by the imaging y-ray telescope 



SIGMA aboard the GRANAT spacecraft. These events were in- 
terpreted as arising from positron annihilation near a galactic 
black hole ( [Ramaty & Meszaros|1981); however the IE 1740.7- 
2942 1-day emission in September 1992 was not confirmed by 



OSSE ( Jung et al. 1995). We note that the different physical char- 



acteristics of the environment inferred from the line width mean 
there is no simple connection between these observations and 
the diffuse 511 keV galactic emission. 

Several scenarios have been proposed to explain observa- 
tions. In general, because of the poor angular resolution of SPl, 
we cannot exclude the possibility that the Galactic 511 keV 
gamma-ray line emission originates in a limited number of point 
sources ( |Knodlseder et al. 2005| l. However, the spectral prop- 
erties of the 511 keV line measured by SPI ( Jean et ar]|2004| l 
suggest that a large fraction of positrons generated in compact 
objects should propagate in the interstellar medium before being 
annihilated. The SPl data can be explained in terms of the anni- 
hilation of positrons in a static warm (~ 10'* K) or cold (~ 10^ K) 
slightly (a few %) ionized medium or, as has been recently pro- 
posed by Churazov et al. (2010 1, with a model of positron anni- 
hilation in a radioactively cooling interstellar medium. Among 
the compact objects, we know that the supernovae must emit 
positrons by the decay of the nuclear synthesis products. One 
problem in this case is that the bulk of positrons are observed 
from the Galactic bulge, in contrast to what we would expect 
from the supernovae distribution in the Galactic plane. A pos- 
sible solution has been proposed by Higdon et al. (20091, who 
suppose that a large amount of positrons generated in the Disk by 
the supernovae explosions are transported to the bulge following 
the galactic magnetic field lines. The low mass X-ray binaries 
(LMXBs) have been proposed by |Weidenspointner et al. ( |2008a| l 
as a class of positron sources that could explain the 511 keV 
emission line measured by INTEGRAL. Basically, this idea is 
supported by two arguments. First, the LMXBs are old systems 
concentrated in the Galactic bulge, where we detect the main 
amount of annihilation radiation. The second reason is that the 
SPI results show an asymmetric 511 keV emission morphology 
with a spatial shape that is consistent with the LMXBs distribu- 
tion above 20 keV. In X-ray binaries, positrons in e^e^ pairs can 
be created in the vicinity of the compact object, either in the hot 
inner accretion disk, in the X-ray corona surrounding the disk, 
or at the base of the jet, which may channel a fraction of e^e^ 
pairs out of the system (Prantzo s et al.|20T 0'). Finally, the anni- 
hilation can be also explained by new astrophysical ideas. The 
possibility that the annihilation radiation is due to the e^e" pairs 
produced when X-rays generated from the Galactic center black 



hole interact with ~ 10 MeV temperature blackbody low-mass 
(10'^ g) black holes within 10'"^ - lO'^ cm of the center is dis- 
cussed by Titarchuk & Chardonnet ( 2006[ ). Annihilation of light 
dark matter (DM) particles into e*e~ is also a possible explana- 
tion ( IBoehm et al.,2004) ). 

In this paper we use IBIS (the y-ray Imager on Board the 
INTEGRAL Satellite) to estimate the 511 keV flux upper limits 
on possible Galactic point sources. In ^we describe the data 
analysis, and in ^and ^we present and discuss our results. 



2. Data analysis 

2.1. Data sample 

The last session of our IBIS 511 keV data reduction started in 
April 2008. Our data set consists of all the IBIS public data 
available at that time, i.e. about 5 years of observations, from 
October 17, 2002, when INTEGRAL was launched, until April 

2007. In addition, we used all the Core program data until April 

2008. The data sample, which corresponds to the 4th IBIS hard 
X sky survey ( [Bird et al.|201 0') has been filtered to exclude the 
observations affected by strong background and systematic noise 
(typically due to solar ffares). 

All the data that we reduced correspond to 39413 IBIS point- 
ings, also called science windows (ScWs), each lasting about 
2000 s. The total observing time in our data sample is ~ 80 Ms. 
The sky exposure map we obtain with this data set is shown in 
Figure [T] We note that the maximum exposure (about 10 Ms) 
corresponds to the Galactic center region, the region where the 
bulk of the positron emission has been detected by the SPI spec- 
trometer aboard flie INTEGRAL satellite. 



2.2. Standard analysis 

The IBIS telescope is based on the use of the coded aperture 
mask technique ( [Caroli et al. 1987 1, that is a key method for 
providing images at energies above tens of keV, where photon 
focusing becomes impossible using standard grazing techniques. 
This established technique is remarkably successful for weak ex- 
tragalactic sources, as well as for crowded galactic regions like 
the Galactic center In a coded mask telescope, the incoming 
source photons are, before detection, encoded by a mask, so the 
image of the sky has to be reconstructed by decoding the ob- 
servation afterward. The coded masks optimize the background 
subtraction capability because of the opportunity to observe the 
source and the off-source sky at the same time. In fact, for any 
particular source direction, the detector pixels are divided into 
two sets: those capable of viewing the source and those for which 
the flux is blocked by opaque mask elements. The coded mask 
telescope are also characterized by a very large field of view 
(FOV). The FOV is divided in two parts, the fully coded (EC) 
FOV for which all source radiation directed towards the detec- 
tor plane is modulated by the mask and the partially coded (PC) 
FOV for which only a fraction of it is modulated by the mask. 

The IBIS data and software are delivered by Qie INTEGRAL 
Scientific Data Center (ISDC) in Geneva ( |Courvoisier et alT] 
2003). The Off-line Scientific Analysis (OSA) software is reg- 
ularly released by ISDC with bug fixes and some improvement, 
we used the release 7.0 of this software in our data analysis. 
The scientific data analysis methods, hereafter called Standard 
Analysis, are reported by.Goldwurm et al. (2003). In this work 



the IBIS data have been reduced using the lASF-INAF comput- 
ing facility in Rome ( [Federici et al.|20(j9] l. 
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2.3. The IBIS Monte Carlo 

The first motivation behind the effort by the community of IBIS 
in developing the Monte Carlo is to calculate the detector re- 
sponse matrix. The response matrix is a table that determines 
the probability p(E, k) that a photon of energy E, incident on the 
detector plane, produces a signal in the detector channel k. The 
IBIS modeling activities (Laure iiFet al.|[2003,) were performed 
step by step by firstly constructing a calibrated model of each 
IBIS sub-system: ISGRI, PICsIT, Mask, active shield, etc. A 
calibrated model is defined as a GEANT Monte Carlo model, 
further checked by intensive corresponding calibration. The cal- 
ibrated models have been integrated in a second step in the whole 
IBIS mass model. This mass model was then checked using the 
data obtained during the IBIS telescope calibration (Bird et al. 
|200 3 ) in the payload ground calibration at ESA/ESTEC, and re- 
fined using in-flight calibration data (Terrier et al. 2003 | l acquired 
during the whole INTEGRAL mission. 

In this work, we apply the Monte Carlo to the evaluation of 
the detector effective area at 5 11 keV, in order to convert the 
detector counts in photons per unit of area. 

2.4. Evaluation of the 51 1 keV line flux upper limit 

The y-ray sensitivity is defined as the minimum flux detectable 
at a level that can be distinguished from the noise; the lower the 
sensitivity, the higher the signal-to-noise ratio for a given source. 
The sensitivity depends on the significance level, so we need to 
define at what level it becomes a reliable signal coming from a 
real source. In this work we assume that we get a positive signal 
if it is measured at a significance level greater than 2 cr. 

Irrespective of the imaging system we use, the detector count 
rate C is the sum of the source S and the instrumental back- 
ground B. Then 



C-B 



(1) 



For any given source, in a coded mask system, the source and 
background counts can be measured using the detector areas that 
coiTespond to open and closed mask elements, respectively. In 
terms of count rate, the astrophysical measures in the hard X and 
y-ray domains are dominated by the background ( B » 5). In the 
Standard Analysis, assuming a Poissoniam noise, the variance in 
each reconstructed sky pixel of the FCFOV is given by the total 
count recordered by the detector ( jGoldwurm et al.|2003| l. Thus, 



if a source produces the count rate S on the detector, for a given 
exposure T, the signal to noise ratio is simply 



S XT 



s Source Count 



n ^jTotal Count ^J(S +B)xT 



S X 



(2) 



The source signal S depends on the detector area and on 
the imaging efficiency; i.e., for a given source, the fraction of 
detected photons that are effectively used to form the signal S . 
In general, due to finite detector space resolution and the mask 
transparency, this number is lower than one. More precisely, 
the imaging efficiency depends on the ratio between the mask 
pixel and the detection units size (the more pixels, the higher 
the imaging efficiency) and on the real transparency of the open 
and closed mask elements. The mask thickness of 16 mm guar- 
antees a very low transparency, to -3% of the close elements 
at 511 keV, and a high transparency, t[ = 90% of the open el- 
ements. The measured values of the open-mask element trans- 
parency versus the incident angle are reported for many source 
energies by [Sanchez et al.| ( |2005| l: at energies higher than ~ 
30 keV, in particular at 5 1 1 keV, the transparency has a low de- 
pendence on the incident angle. The imaging efficiency effimg 
can be evaluated by the expression ( |Skinner|2008| l 

U(l-i)(fi-fo)V4/(l-/) if^ 

where the mask open fraction / is 0.5 for IBIS, d and m are the 
size of the detector and mask units. For IBIS/ISGRI, because 
d/m - 4.6/11.2 - 0.4107, taking the mask's closed and open 
element transparencies to and ti into account, the imaging effi- 
ciency is 



< m 



> m 



(3) 



effi,„g(5nkeV)^0.75 



(4) 



3. Results 

3.1. ISGRI background and effective area at 511 keV 

As discussed in the previous section, the IBIS sensitivity depends 
on the instrumental background and the detector's effective area 
at 5 1 1 keV. Since the background is caused by radiation induced 
by cosmic rays interactions with the satellite, the telescope, and 
the upper atmosphere, the background is anti-correlated with the 
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Fig. 2 The 511 keV line background as measured by ISGRI. The energy band width is 1 FWHM of the detector energy resolution at 
511 keV. Each INTEGRAL revolution in the horizontal axis lasts 3 days. 
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Fig. 3 A simulation of a 51 1 keV point source as seen by IBIS/ISGRI during a pointing in the early (rev. 50), left panel, and late (rev. 
600), right panel, orbits of the INTEGRAL satellite. The signal level of the detected 511 ke V source goes to 10 cr in the first image 
to 7 cr in the second one. These images have been obtained adding to the real data the data generated by the IBIS Monte Carlo and 
then processing the result with the standard analysis software until the image level. 



solar activity, reaching a minimum at the solar max in the early 
periods of the INTEGRAL activity. During this period, we mea- 
sured a 511 keV backgrouncQ of about % counts Is. The long- 
term background trend is characterized by a general increase 
with time and a mean value of about 10 counts /s (Figure |2|i. 
The peaks in this curve are due to the high-energy component 
of solar flares, which produce bursts of radiation across the elec- 
tromagnetic spectrum from radio waves to X and y-rays. In our 
data sample, we removed any time intervals associated to these 



' We define the 511 keV background as the count rate in a FWHM 
width bin centered at 5 1 1 keV 



bursting events, because they are affected by systematic noise, 
so are not useful for scientific data analysis. 

The ISGRI effective area is estimated by the Monte Carlo 
method (5 2.3i. A mono-energetic 511 keV point source is ob- 



tained by simulating a parallel beam incident on the top of the 
detector. If A^„, is the number of input photons, the effective area 
is simply estimated by 



A,f f (511 keV) ^ Abeam y- 



N, 



det 



(54.2 + 0.2) cm^ 



(5) 



where Aheam is the area covered by the beam, and Njei is the 
number of the photons detected in a FWHM width bin centered 
at 511 keV. Taking the ISGRI geometrical area (~ 3000 cm^) 
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Crab 

a 




Fig. 4 The Crab nebula as detected by IBIS/ISGRI in the soft y-ray domain with a 3 Ms of exposure. Left and right top panel: the 
ISGRI 40-100 keV shadowgram accumulated during one ScW (no. 010200210010) and the decoded sky image. Left bottom panel: 
sky mosaic in 431-471 keV (colorbar: -2-8 cr). Right bottom panel: sky mosaic in 491-531 keV (colorbar: 3 - 4cr). 



Table 1 The Crab nebula as observed by IBIS with an exposure 
of 3 Ms. 



Energy (keV) 


signal (tr) 


counts/s 


20 - 40 keV 


8372 


157.527 ±0.019 


40- 100 keV 


5360 


93.984 ±0.017 


100 - 300 keV 


1250 


17.608 ±0.014 


431 -471 keV 


9.5 


(3.0 ±0.3) 10-2 


491 -531 keV 


3.8 


(1.2 ±0.3) 10-2 


551 - 591 keV 


not detected 





into account, the effective area results from of an about 2 % peak 
efficiency at 5 1 1 ke V. 

3.2. IBIS 511 keV line sensitivity 

If the 511 keV flux measured by SPI was due to a set of point 
sources, would it be possible to detect some of these sources 
with IBIS? To answer this question we consider the simple case 
where all the flux is concentrated in just one point source. In this 



parameter 


value 


photon index for E < Etrmk 
Ebn-ak (keV) 

photon index for E > Etrmk 
photons crrT^ .s"' keV^^ at 1 keV 


2.08 ± 0.01 
100 (frozen) 
2.23 ± 0.03 
9.3 ± 0.2 



Table 2 The broken power law model for the Crab spectrum in 
X- and soft y-rays reported by ,Jourdain et al. (2008) 

case, we would detect a flux of phcmT^ i"' with a count 
rate of 

Siuk.v = 10-3 X ^ X "^ffimg = (2.03 ± 0.01) X 10"^ s'' (6) 

From equation (j2]i it follows that for an exposure of T = 10^ s, 
assuming a background Bjiifoy — IQcounts/s, we obtain a sig- 
nal to noise of 

- =55iifeyX J-^— =20cr (7) 

n y 05l\keV 
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Table 3 511 keV 2 a Sgr A* upper flux limit in the Galactic Center visibility periods. 



Period 


Start Time 
(UT) 


End Time 
(UT) 


INU. vv s 


ivis 


nUA IlllllL 

lu pnciii i 


opung 


2003 


-02 


28 


2003 


04-23 


1 / J 1 


U.OJU 


J. I 


2003 Autumn 


2003 


-08 


10 


2003 


10-14 


1717 


1.872 


3.2 


2004 Spring 


2004 


-02 


16 


2004 


04-20 


1862 


1.046 


4.5 


2004 Autumn 


2004 


-08 


17 


2004 


10-27 


2191 


1.292 


4.3 


2005 Spring 


2005 


-02 


16 


2005 


04-28 


2144 


1.165 


4.7 


2005 Autumn 


2005 


-08 


16 


2005 


10-26 


1667 


0.790 


5.9 


2006 Spring 


2006 


-02 


09 


2006 


04-25 


1869 


1.501 


4.1 


2006 Autumn 


2006 


-08 


16 


2006 


11-02 


1770 


1.080 


5.0 


2007 Spring 


2007 


-02 


01 


2007 


04-22 


985 


0.347 


9.2 



Notes. 

Each ScW lasts about half an hour. 

**' Obtained from the exposure map. The exposure map, and in particular the value in the Galaxy depends both on the number of ScWs and the 
dithering patterns 



Table 4 511 keV 2 cr flux upper limits for the hard X (E > 20 keV) microquasars reported in the 3rd IBIS catalogue (Bird et al. 
[2007] l. 



Source name 


R.A. 


Dec 


Error" 


-r 20-40 


-f^ 40- 100 


2 a flux limit 


Type'' 




(deg) 


(deg) 


(arcminutes) 


(mCrab) 


mCrab 


(10-* phcm-^ J-') 




XTE J1550-564 


237.745 


-56.479 


0.2 


34.0 ±0.1 


55.3 ±0.2 


2.7 


LMXB 


Sco X-1 


244.980 


-15.643 


0.2 


685.7 ±0.3 


24.7 ± 0.3 


3.7 


LMXB 


GRO J1655-40 


253.504 


-39.846 


0.6 


2.3 ±0.1 


2.7 ±0.2 


2.5 


LMXB 


GX 339-4 


255.706 


-48.792 


0.3 


40.7 ±0.1 


46.7 ± 0.2 


2.6 


LMXB 


IE 1740.7-2942 


265.978 


-29.750 


0.2 


29.8 ±0.1 


36.6 ±0.1 


1.6 


LMXB 


GRS 1758-258 


270.303 


-25.746 


0.3 


58.8 ±0.1 


75.3 ±0.1 


1.6 


LMXB 


SS 433 


287.956 


4.983 


0.5 


10.4 ±0.1 


5.2 ± 0.2 


2.4 


HMXB 


GRS 1915+105 


288.799 


10.944 


0.2 


296.8 ±0.1 


123.4 ±0.2 


2.5 


LMXB 


Cyg X-1 


299.590 


35.199 


0.2 


763.7 ±0.2 


876.7 ±0.3 


3.0 


HMXB 


Cyg X-3 


308.108 


40.956 


0.2 


196.5 ± 0.2 


78.3 ±0.3 


2.6 


HMXB 



Notes. 

Position error expressed as radius of 90 % confidence circle in arcminutes. 

Fluxes averaged on the total exposure: 
20-40 keV: 1 mCrah = 7.57 x 10^'^ erg cm-^ .s"' = 1 .7 1 x lO^* ph an^ s'^ 
40-100 keV: 1 mCrab = 9.42 x lO"'- ergcrrr^ i-' = 9.67 x 10"^ phcm.-^ i"' 

Type identifiers: LMXB=Low Mass X-ray binary, HMXB=High Mass X-ray binary. 



Table 5 511 keV 2 cr flux upper Umits for the hard X (E > 20 keV) Supernova Remnants (SNR) reported in the 3rd IBIS catalogue 
dBird et al.|2007] l. 



Source name 


R.A. 


Dec 


Error" 


-^20- 


b 

-40 


-^40- 


b 

100 


2 a flux limit 




(deg) 


(deg) 


(arcminutes) 


(mCrab) 


mCrab 


(10^^ phcm-\?-') 


4U 0022-f63 


6.319 


64.159 


3.6 


0.7 ± 


0.1 


0.7 ± 


0.2 


2.3 


PSR J1811-1926 


272.827 


19.417 


2.9 


0.8 ± 


0.1 


1.0 ± 


0.2 


14 


IGR J18135-1751 


273.395 


-17.871 


2.2 


1.1 ± 


0.1 


1.7 ± 


0.2 


2.3 


SNR 021.5-00.9 


278.388 


-10.579 


1.2 


3.3 ± 


0.1 


3.3 ± 


0.2 


3.0 


AX J1838.0-0655 


279.509 


-6.916 


1.5 


1.9 ± 


0.1 


3.0 ± 


0.2 


2.9 


Kes 73 


280.338 


-4.948 


1.3 


2.2 ± 


0.1 


4.2 ± 


0.2 


2.9 


AX J1846.4-0258 


281.596 


-2.983 


1.8 


1.7 ± 


0.1 


2.4 ± 


0.2 


2.8 


Cas A 


350.848 


58.815 


0.9 


4.1 ± 


0.1 


2.4 ± 


0.2 


2.5 



Notes. 

Position error expressed as radius of 90 % confidence circle in arcminutes. 

Fluxes averaged on the total exposure: 
20-40 keV: I mCrab = 1 .51 x 10"'^ ergcra^^r' = 1.71 x IQ-* phcmT^ s'^ 
40-100 keV: 1 mCrab = 9.42 x 10"'" ergcm-^ i"' = 9.67 x 10"' phcm-^ .s"' 



We conclude that for a 10 Ms exposure the 2cr 51 1 keV sensitiv- 
ity is on the order of lO"'* ph cm^^ s^\ 

The IBIS sensitivity can also be estimated by using data 
from the Crab observations. The Crab Nebula was monitored by 
INTEGRAL through a large number of pointings that were per- 



formed during many dedicated observations. The Crab exposure 
that we obtained with these data is equal to 3 Ms. Reducing all 
the available data, we detect in the sky mosaic the Crab Nebula 
from 20 keV up to 531 keV (Fig.|4]i. The signal goes from 8372 cr 
in the 20 - 40 keV band to 3.8 cr in the 491 - 531 keV (Table 
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LMXBs 511 keV flux upper limits 



S 14 



h2 

10 



■ 



I 



nr. 



1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 5.5 6.0 6.5 7.0 7.5 8.0 8.5 
51 1 keV upper limit (1 OM photons / cm»2 s) 



we call spring and autumn as they occur mainly in these sea- 
sons in the northern hemisphere. A Ms time scale is explored by 
considering the observation of the Galactic center during these 
visibility periods. To select the data, we extracted from the data 
archive all the IBIS pointings with Sgr A* in the field of view. 
We then obtain a 5 1 1 keV sky mosaic for each Galactic center 
visibility period. No significant excess in the 511 keV bin are 
detected. We report the 511 keV flux limits in the center of the 
Galaxy during the INTEGRAL visibility periods in Table [3]. 

3.4. 51 1 keV flux upper limit for compact sources 

The LMXBs detected by IBIS in the Galactic plane could make 
a substantial contribution to the galactic 511 keV emission (see 
Weidenspointner et al.||2008a|b| l. Among a total of 79 LMXBs 



reported in the hard X domain by Bird et al. ( 2007 1, 70 are lo- 
cated near to the Galactic plane (\b\ < 10), and then can in prin- 
ciple contribute to the Galactic 511 keV emission. For each of 
these sources, we evaluate the 511 keV flux limits. Depending 
on the exposure, the flux limit goes from 1.6 x 10^* pit cm^^ i"' 
in the maximum exposure (10 Msec) sky regions to about 4.0 X 
limit of 1.6 X IQ-'^ phcnT^ s-\ estimated for the IBIS sources phcmT^ i^' (Fig.lsl. The only exception is 4U 0614+091, 



Fig. 5 The distribution of 5 1 1 keV flux upper limits for the IBIS 
LMXBs population. The first peak coiTespond to a 2cr upper 



located in the Galactic Bulge 

As the last energy band is where we are looking for a direct 
positron annihilation signal from point sources in the data, this 
signal gives us the opportunity to estimate the IBIS sensitivity at 
511 keV using real data. A recent Crab measurement made with 
the INTEGRAL/SPI y-ray spectrometer is reported by IJourdain 
|et al.| ( [2008] l, where the data are modeled by a broken power law 
(Table^grrhe 491-531 keV Crab flux is equal to 



Fcrab(49l - 531 keV) = (6.8 ± 0.5) x IQ-* p hem 



-2 



(8) 



Since we detect this flux at a significance level of 3.8 cr, the sen- 
sitivity at 2 cr level is 3.6x lO"'* ph cmT^ s"', obtained with an ex- 
posure of 3 Ms. With a 10 Mi exposure the sensitivity (Equation 
^ will be 



S^'fiSllkeV) = (2.0 ± 0.7) x 10">/z cm"^ s 



-2 „-l 



3.3. The Galactic center at 51 1 keV 



(9) 



for which we obtain a higher flux limit (8.7x1 0" ph cmT ), 
because it is located in a sky region with a lower exposure. 

Among X-ray binaries (XRBs), microquasars (MQs) are 
promising galactic positron emitters (see Guessoum et al._2006| . 
We therefore report the 511 keV upper limits for the X-ray 
Binaries classified as microquasars and detected by IBIS above 
20 keV (Table |4|. Among the 10 MQs in our sample, 7 are 
LMXBs, 3 High Mass X ray Binaries (HMXBs). 

Some of the atomic nuclei produced in supernovae produce 
positrons by radioactive fi^ decay of the nucleosynthesis prod- 
ucts. The initial positron energy is expected ~ 1 MeV, in agree- 
ment with the injection energy constraint (|Beacom & Yiiksel 
2006). The 511 keV flux upper limits obtained with our data 
analysis for the supernova remnants (SNRs) detected by IBIS 
above 20 keV are reported in Table[5] These limits have been es- 
timated by considering these objects as possible 511 keV point- 
like sources, i.e. with an extension lower than the IBIS angular 
resolution (12 arcminutes). However, the y-ray emission from 
SNRs could be more extended, in which case the sensitivity 



With the only exception given by the Crab, detected at about 4 
sigma of significance (Table[T]i, with about five years of data, we 
obtain a 51 1 keV (491-531 keV bin) all sky map of our Galaxy 
without finding any signal. We note that, as the Crab signal is 
consistent with the expected continuum contribution at this en- 
ergy, we do not find any evidence of a 5 1 1 ke V emission line in 
this source. At each point on this map the variance is given by 
the Standard Analysis of the raw data. Then, at a given position, 
dividing the standard deviation into the detector effective area 
and imaging efficiency, we estimate the sensitivity at 511 keV. 
In the Galactic center (Sgr A*) position, we estimate for an ex- 
posure of 10 Ms, the following flux upper limit (given at at 2 cr 
of significance level) 



S2^{SgrA*) 



phcnT^ s-^ 



(10) 



We obtain a similar 511 ke V flux upper limit for the sources 
located in the Galactic bulge, where the observations have an 
exposure of about 10 Ms. We then focused our interest on the 
this region on different timescales. Due to the requirement on 
the orientation of the INTEGRAL solar panels, the region of the 
Galactic center can be observed in only two periods per year that 



degrades almost linearly with the source extent (Renaud et al. 
|20U6l ). 

4. Discussion 

We have focused data analysis on the search for 511 ke V line 
emission from point sources. A possible detection at this energy 
would be a proof of an e^e^ annihilation process in Galactic 
sources, providing an important clue to the mystery on the ori- 
gin of Galactic bulge positrons. Through a data analysis of about 
five years of observations, we do not detect any 511 keV y-ray 
point source with IBIS. The flux upper limits that we report for 
the sources located in the Galactic bulge are on the order of 
2 X 10""* ph cmT^ 5"' . A direct consequence of this result is that, 
if the flux of positrons measured by SPI from the Galactic bulge 
(~ lO"-' phcm^^ .?"') was caused just by point sources, assum- 
ing the same intensity for all of them, there should be at least six 
sources. As the sensitivity degrades with energy, a better con- 
straint on the 511 keV emission should be achieved by focus- 
ing data analysis to the lower energy continuum emission due 
to ortho-positronium. However, although the value of fraction 
of positronium measured by SPI for the diffuse Galactic bulge 
emission is close to one, the annihilation radiation from Galactic 
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objects may bias the measurements. In this regard, we note that 
the SPI data analysis is optimized for narrow gamma lines with 
poor sensitivity to broad emission lines. In terms of data anal- 
ysis, the positronium continuum needs to be separated by the 
intrinsic source continuum emission up to a few hundred keV. 
Therefore, to avoid contamination, a simple imaging technique 
can use an energy bin with just a fraction of the positronium 
continuum due to positronium, thereby degrading the sensitivity 
that should be obtained with this method. A different approach, 
based on spectral analysis, may be to search for any positronium 
components in the continuum emission of the Galactic sources. 

Taking the Galactic center distance into account, the mea- 
sured 511 keV flux corresponds to the annihilation of 2 x lO'*^ 
positrons per second. We note that, since the t otal X-ray lu- 
minos ity estimated in LMXBs is ~ 10^^ erg/ s ("Grimm et al.' 
2002 1, the energy (2 x 10^^ ergjs) required to create the Galactic 
positrons can in principle come from these sources. The idea 
that these sources could make an important contribution to the 
annihilation e^e" radiation is also supported by the correla- 
tion between their spatial distribution and the distribution of 



the 511 keV radiation measured by SPI ( Weidenspointner et al. 



2008a| l. The fact that we do not detect any LMXB witii IBIS is 
consistent with the idea that a significant part of the positrons 
propagates in the interstellar medium before being annihilated 
away from the birth place. For the microquasars, a possible 
511 keV annihilation radiation can arise near the source, from 
the interaction of the jet with the interstellar medium, or with 



the companion star in the case of a misaligned jet (Guessoum 
et al.|2006| . In general, a possible electron-positron annihilation 



emission due to a Galactic jet would be seen as a point source if 
its angular size is lower than 12 arcminutes (i.e. the IBIS angu- 
lar resolution). The study of the jets from microquasars makes 
use of radio and X-ray observations (for a recent review see 
Gallo 2010| l. The jet interpretation of the radio emission was 
confirmed by the miUiarcseconds scale jets observed in Cygnus 
X-1 ^Stirling et al.| [200TT) and GRS 1915+105 ( |Dhawan et al"! 
200011. hi X -ray band, arcsecond scale jets have been observed 
for the first time in SS 433 with the Chandra X-ray observatory 
( [Marshall et al.|2002| ). For microquasai's 1 El 740.7-294 2, abroad 



emission line detected by SIGMA has been reported ( Sunyaev 
et al.|199T i. However, since it was not confirmed by OSSE data 
Pung et al.||1995[ ), this measure is quite controversial. With a 



5-year sky monitoring at 511 keV with IBIS, we do not detect 
this phenomenon, showing that, if true, it must have a rather low 
cycle of activity. 

The 511 keV Galactic luminosity could be explained by 
positrons generated by y6^ decay of the nucleosynthesis products 
in supemovae. If a significant fraction of positrons are annihi- 

( 2010) have proposed that 



lated near the birth place, Martin et al. 



supernovae could be detected as 51 1 keV point sources, but the 
flux is expected to be lower than the actual instrument sensitivity 
at this energy. Depending on the temperature of the environment 
where the annihilation takes place, a 511 keV broad emission 
line cannot be ruled out for supernovae. If the process takes place 
in the hot medium like the shocked regions of Supernova rem- 
nants, the expected line broadening would be ~ 30 keV for a 
temperature of 10^ K, and then the emission would be within the 
energy bin (491-531 keV) used in our data analysis. 
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